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CN ■ Abstract 

Dynamics and stability of average current control of DC-DC converters are analyzed by sampled- 
| data modeling. Orbital stability is studied and it is found unrelated to the ripple size of the 

orbit. Compared with the averaged modeling, the sampled-data modeling is more accurate and 
systematic. An unstable range of compensator pole is found by simulations, and is predicted by 
sampled-data modeling and harmonic balance modeling. 
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1 Introduction 



In this paper, dynamics and stability of average current control (ACC) of DC-DC convert- 
q | ers are analyzed by sampled-data modeling and harmonic balance modeling. Previous works 

[U [21 O HI [6] on ACC DC-DC converters are generally based on averaged modeling. In the 
averaged modeling, each circuit module is modeled separately and approximately. For example, 
the switches are approximated as a three-terminal model. The current loop is obtained using 
sampled-data dynamics and then approximated by continuous-time dynamics. The PWM mod- 
ulator is approximated as a constant gain. The duty cycle, a discrete-time variable, is treated 
as a continuous-time variable. Then all of these approximate modules form a continuous-time 
model. Therefore, the orbital nature of the nominal periodic solution is lost. Instead, an equi- 
librium is obtained as the nominal solution, and graphical ripple analysis is used to determine 
stability. 

In contrast, no such approximations are involved in the sampled-data modeling. The or- 
bital nature of the nominal periodic solution is preserved, and exact orbital stability can be 
determined. Also, the sampled-data modeling focuses on system operations, especially switching 
actions. The switching action is very important for the derivation of system dynamics, and it is 
accurately formulated in the sampled-data modeling. Therefore, the advantage of the sampled- 
data modeling is that it is more accurate and systematic. Although sampled-data analysis of 
converters has been a topic of investigation for the past two decades El El [101 [TTJ , this pow- 
erful tool is not widely used. This paper has an aim to increase the appreciation and use of the 
sampled-data modeling. 

Harmonic balance is another approach to analyze subharmonic oscillation (period-doubling 
bifurcation) in DC-DC converters [12]. By using this approach, a simple analytical form for 
stability is obtained and it implies existence of an unstable range of compensator pole, which is 
confirmed by simulations. 

The remainder of the paper is organized as follows. In Section [21 the operation of average 
current control is reviewed and an instability example is given. In Section [31 a general sampled- 



data model for DC-DC converters developed in [13] is reviewed and applied to the ACC DC-DC 
converter. In Section U harmonic balance analysis developed in |12j is applied to study period- 
doubling bifurcation in the ACC buck converter. Conclusions are collected in Section 



2 Operation and Instability of Average Current Control 

An ACC DC-DC converter is shown in Fig. [TJ Let switching frequency f s = 1/T, and ui s = 2irf s . 
The ACC operation is as follows [3]: The inductor current ii is sensed by a resistor R s and 
compared with a current reference v r from a voltage loop (not shown). The difference is amplified 
by a compensator. A typical compensator is 
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The compensator output y(t) is compared to a T-periodic ramp signal h(t), which has h(0) = V{ 
and h(T~) = Vh, to generate a switching signal to the power stage. The nominal steady-state 
solution of a DC-DC converter is a T-periodic orbit, not an equilibrium point as depicted in the 
averaged modeling. Denote the orbit by x (t), and the corresponding steady-state compensator 
output by y°(t). An orbit (like x°(t) or y°(t)) is (orbitally) stable if a state trajectory stays in 
the orbit if it is perturbed. If the state trajectory eventually leaves the orbit, then the orbit is 
unstable. Previous works [H El IH E] imply that u p and ripple size of y affect stability linearly. 
From the following simulations, however, an unstable range of u p is found and the ripple size is 
unrelated to orbital stability. 




Figure 1: System diagram of an ACC DC-DC converter 

Example 1 Consider an ACC buck converter [2J p. 114]. The system parameters are v s = 14 V, 
nominal (set-point) output voltage Vset = 5 V, v r = 0.5, f s = 50 kHz, L = 46.1 juH, C = 380 
(jF (with ESR R c = 0.02 fi), R = 1 Q, R s = 0.1 Q,Vi = 0, V h =l, K c = 75506 and u z = 5652.9 
rad/s. 

Compensator pole oj p is varied from 0.14w s to 0.81cj s . The compensator output ripple size 
increases as uj p increases. For lo p < 0.19w s , x°(t) or y (t) is stable. For 0.19o; s < co p < 0.490^, 
x°(t) is unstable. Take u p = 0.21a; s for example, y°(t) is unstable with small ripple (Fig. [2]). For 
uj p = 0.49w s , period-doubling bifurcation occurs, and the unstable T-periodic orbit (Fig. [3|) and 
the stable 2T-periodic orbit (Fig. U|) coexits. The two orbits have the same averaged duty cycle 
0.357 and share the same averaged orbits. Based on the averaged model, they should have the 
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same stability. However, their stabilities are different. For uj p > 0A9u s , x°(t) or y°(t) is stable 
again. Take u p = 0.81o; s for example, y°(t) is stable with large ripple (Fig. [5]). □ 
The orbital stability and the unstable range of oj p will be explained by the sampled-data 
modeling and harmonic balance modeling discussed next. 
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Figure 2: Unstable y°(t) with small ripple, u p = 0.21u; s 
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Figure 3: Unstable y°{t), lo p = 0.49w, 



3 Sampled-Data Analysis 

The operation of average current control can be described exactly by a block diagram model 
developed in [13], which studied closed-loop DC-DC converters. Here the model is applied to an 
ACC DC-DC converter, with the voltage loop open and the reference signal v r considered as a 
control variable. 
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Figure 4: Stable 2T-periodic orbit, oj p = 0.49oj s 




Figure 5: Stable y°(t) with large ripple, co p = 0.81a; 
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The block diagram model is shown in Fig. [6l and it can model many DC-DC converters, 
such as buck, boost, and Cuk converters, etc. In the model, A\,A2 E H NxN , B±,B2 E R Arx2 , 
C,E 1 ,E 2 E R lx7V , and D E R lx2 are constant matrices, where N is the state dimension, 
typically given by the number of energy storage elements in both of the power stage and the 
compensator. For example, N = 3 for a typical second order buck converter with a first order 
compensator. Within a clock period, the dynamics is switched between two stages, S% and 
S2. Switching occurs when the ramp signal h(t) intersects with the compensator output y := 
Cx + Du E R. The system is in S\ when y(t) > h(t), and switches to S2 at instants when 
y(t) < h(t). Here the switching operation is modeled exactly, without seperate and graphical 
approximations as in the averaged model. 

h(t) = V l + (V h -V l )(±modl) 
clock 

y = Cx + Du 



Switching 
Decision 



Switch to S*i or £2 



Si : < 
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= E x x 


S 2 : < 
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= A 2 x + B 2 u 




= E 2 x 



Figure 6: Block diagram model for ACC DC-DC converter 



Example 2 Consider Example 1. Let the state x = (*£,, vc, v e i, v e2 )' , where il is the inductor 
current; vc is the capacitor voltage; v e \ and v e % are the states of the compensator. Then, 
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□ 

Let x n = x(nT). Assume the switching frequency is high enough that u = (v s ,v r )' E R 2xl 
can be considered constant within the cycle, and denote its value by u n = (v sn ,v rn )' , which can 
be varied from cycle to cycle. The notation v sn , instead of v s ^ n , is used for brevity. 

The periodic orbit x (t) in Fig. O corresponds to a fixed point x (0) in the sampled-data 
dynamics. Let the steady-state duty cycle be D c and d = D C T. A typical periodic solution x°(t) 
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is shown in Fig. where x {d ) = Aix°(d) + B\u and x°(d + ) = A 2 x°(d) + B 2 u denote the time 
derivative of x°(t) at t = d~ and d + , respectively. 



;,0 



(d+) 




x°(d) 



Figure 7: A typical periodic solution of a DC-DC converter in state space 

Using a hat * to denote small perturbations (e.g., x n = x n — x°(0)). From [13], the linearized 
sampled-data dynamics is 

X n+1 = <M n + Tu n = <M n + TiV sn + T 2 V rn (2) 

where $ e R^^ and T = [T 1 ,T 2 ] G K Nx2 are 

$ = e ^(T-d) a (0*1 ~ A 2 )x°(d) + (gi - B 2 )u)C Ald 

C(A 1 x°(d) + B 1 u) -h{d) 

= ^ 2 (T- d ) fJ {x»{d-)-x\d+))C Md 

K Cx°(d-)-h{d) ' V ; 

T = e A ^ T ~ d H ( d e A ^daB l - ~ ^ } (C f V° daB, + D)) + f"" e A ^daB 2 

Vo Ca; (d-) - V Jo "Jo 

Local orbital stability of the converter is determined by the eigenvalues of The periodic 
solution x°(t) is asymptotically orbitally stable if all of the eigenvalues of $ are inside the unit 
circle of the complex plane. 

Example 3 Consider again Example 1. The eigenvalues of <3? are calculated as co p is varied from 
0.14w s to 0.81UV For uj p = 0.14w s , all eigenvalues are inside the unit circle and x°(t) is stable. 
As Up is increased, the eigenvalue trajectory leaves the unit circle through -1 when u p = 0.19u; s , 
and enters again the unit circle when uj p = 0.49w s . This explains exactly the unstable range of 
ojp reported in Example 1. □ 

Example 4 For comparison, the linearized dynamics of state-space averaged model of Example 
1 is derived: 

x = (A 1+ Vh _^ )x (4) 

For Up varied from 0.14a; s to 0.81w s , the poles of (HJ) are in the left half of complex plane and 
the converter should be stable according to the average model. This contradicts the simulation 
results in Example 1. □ 
The dynamics © can be also used to derive some useful transfer functions. The reference 
signal v r is used as a control variable. Let E = {E\ +E 2 )/2. From Eq. ([2]), the control-to-output 
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Figure 8: Unified dynamic model of buck converter 



voltage transfer function is 



T oc {z) = ^44 = E(zl - fc)" 1 ^ (5) 
v r {z) 



Given a transfer function in z domain, say T(z), its effective frequency response is T(e^ LuT ), 
which is valid in the frequency range \u)\ < ^. 
Also from Eq. (J2j) , audio-susceptibility is 

T os (z) = V -^\ = E(zI-<S>y 1 T l (6) 

Other transfer functions, such as control-to-inductor current transfer function and output impedance 
can be derived similarly. 

4 Harmonic Balance Analysis of Period-Doubling Bifurcation 

The unstable range of lo p can also be predicted by harmonic balance modeling. In [12], harmonic 
balance modeling is applied to study period-doubling bifurcation of a buck converter under 
voltage or current mode control. The buck converter model is shown in Fig. [HJ Here, the same 
model can be applied to study an ACCbuck converter, with G(s) = R s H c (s)Gi(s), where Gi(s) 
is duty-cycle-to-inductor-current transfer function 

RCs + 1 



RLCs 2 + Ls + R 



Let the critical source voltage be V* such that period-doubling bifurcation occurs for v s > V* . 
A larger V/ implies a larger stable operating range of source voltage. According to [12], a good 
estimate of V* is 

V* « Yhzli : ( 8 ) 

2Re[G(ju s ) - G{^)] 

For an ACC buck converter with oj s 3> 1/VLC and 1/RC, Eq. ([8]) can be simplified as 

V * W 3R~K C m (9) 

where k = uj p /uj s and 

= < 1+fc2 '(^ + " 2 ) (10) 

k 

The function <p{k) = (p(u} p /u s ) is shown in Fig. 0and has a minimum at oj p = 0.38o; s . From @, 
the following operating conditions are desired to have a larger stable range of source voltage: 
larger <p, ramp amplitude (Vh — Vi), L, uj z , oj s ; or smaller gains R s and K c . This gives insight 
on how these parameters affect (period-doubling) instability. A guideline of previous works [5] 
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on ACC is Up = 0.5cv s . From Fig. [9l this choice of uj p results in smaller <j> and smaller stable 



operating range of source voltage. Setting oj p 
larger operating range. Let 



instead (as an example) will have two times 



VI' 



V. 



s lfc=0.38 



0.79(V fe - Vi)Luj z uj s 
R,K r . 



(11) 



In 0, V* is a function of u p . If there exists a oj p such that VQ* < v s , then period-doubling occurs. 
Therefore, a larger vy™ n is desired and the relationship between V™ m and v s determines whether 
there exists an unstable range of u p : 

Theorem 1 Consider an ACC converter with a source voltage v s . The voltage loop is open and 
the current loop has a compensator with a pole uj p and a zero. If v s > V™ m , then there exists 
an unstable range of co p which results in period- doubling instability. If v s < V" im , then u p does 
not affect the stability. 




Figure 9: Function <f>{k) = 4>(uj p /uj s ) 



These two cases are illustrated in the following two examples. 

Example 5 In Example 1, V™ m = 8.57 < v s , therefore an unstable range of uj p exists which 
results in period-doubling. This agrees with the simulation in Example 1. □ 
From Examples 3-5, both sampled-data and harmonic balance models, but not the averaged 
model, predict accurately about period-doubling. 

Example 6 Consider another ACC buck from O p. 982]. The system parameters are v s = 5 
V, Vset = 2 V, v r = 0.279, f s = 180 kHz, L = 13 /xH, C = 750 //F (with ESR R c = 5 mfi), 
R = 0.43 Q,, R s = 0.06 n, V t = 0, V h =2.7, K c = 98000, and uj z = 6723 rad/s. 

Here y^™" = 35.86 > v s and period- doubling does not occur for any value of uo p . However, 
other kind of instabilities may exist. Generally, u p is greater than oj z . To show rich dynamics 
in the converter, let u p = 5655 < uj z . The nominal y°(t) is shown in Fig. [T0| A pair of 
complex conjugate eigenvalues of the sampled-data model are outside the unit circle, which 
implies occurrence of a Neimark bifurcation [llj . Therefore, the orbit of y (t) is unstable although 
it has very small ripple. Also, a pair of eigenvalues of the state-space averaged model ([2]) are 
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calculated and they are found on the right half plane. Therefore, in this example, the averaged 
model dU) also predicts the Neimark bifurcation correctly. Generally, the Neimark bifurcation, 
but not the period-doubling bifurcation, can be predicted by the averaged model. □ 



Figure 10: Unstable y (t) with almost no ripple 



5 Conclusion 

Sampled-data modeling is applied to DC-DC converters under average current control. The 
orbital nature of the nominal periodic solution is preserved. Orbital stability is studied and 
is unrelated to the ripple size of the orbit. Various transfer functions are derived. In the 
sampled-data model, the system dynamics is derived directly and exactly from discrete switching 
operations. In the averaged model, however, each circuit module (e.g., the switch, PWM modu- 
lator and current sensing) is modeled separately and approximately. Therefore, the sampled-data 
modeling is more accurate and systematic than the averaged modeling. 

An unstable range of compensator pole which results in period-doubling instability is found 
by simulations, which can be predicted by the sampled-data modeling and harmonic balance 
modeling. This nonlinear effect can not be predicted by the averaged modeling. How different 
system parameters affect stability is derived. For example, a critical source voltage V* is derived. 
The relationship of its minimum value VJ nm with the source voltage determines whether an 
unstable range of compensator pole exists (Theorom 1). 

Two typical instabilities in DC-DC converters are period-doubling bifurcation and Neimark 
bifurcation, and both can be accurately predicted by the sampled-data modeling. The period- 
doubling bifurcation can be also predicted by the harmonic balance modeling. The instability 
predicted by the averaged modeling is generally the Neimark bifurcation, but not the period- 
doubling bifurcation. 
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